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A problea  of  coating  earth  resource  survey 
aircraft,  proposed  by  NASA,  is  focaulated  as  a 
traveling  salesaan  problea,  in  vhich  the  salesaan 
(aircraft)  has  a range  constraint.  A heuristic 
algoritha  is  presented,  which  seeks  a ainiaal  length 
set  of  subtours  through  n cities.  The  aircraft  begins 
a subtouc  at  the  base  location,  visits  a subset  of  the 
n cities  and  returns  when  the  range  constraint 
prevents  a visit  to  another  city.  Additional  subtours 
ace  created  until  all  cities  are  visited. 


The  algoritha  is  pcogcaaaed  in  FOBTBAN  for  use  on 
digital  conputers.  The  IBH  360/67  coaputec  at  the 
Naval  Postgraduate  School  was  used  to  find  solutions 
to  three  operational  probleas  of  size  seven,  eighteen 
and  twenty-five  cities.  Coaputation  tines  for  each 
problen  was  under  20  seconds  and  the  solutions  were 
significiantly  better  than  feasible  solutions 
calculated  without  the  use  of  the  algoritha. 
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I.  IMTHODOCTION 


The  Natioaal  Aeronautics  and  Space  Administration  is 
fostering  many  technology  utilization  programs  to  explore 
the  application  of  its  aerospace  technology'  to  other  fields. 
Suggesting  how  our  finite  natural  resources  can  be  used  most 
effectively  is  one  of  its  more  important  considerations. 
NASA  has  been  working  to  sophisticate  the  remote  sensing  of 
the  earth  from  a variety  of  aerospace  platforms.  The 
information  available  from  such  a program  can  provide 
accurate  and  timely  data  to  assist  in  the  management  of  the 
earth's  natural  and  cultivated  resources.  NASA  has 
demonstrated  the  usefulness  of  remote  sensing  to  a wide 
variety  of  users  with  emphasis  on  the  disciplines  of 
geology,  agricnlture,  forestry,  oceanology,  meteorology  and 
hydrology.  For  example,  an  agriculturist  can  be  provided 
with  identification  and  area  measurements  of  major  crop 
types  for  use  in  forecasting  potential  crop  yield.  In 
addition,  remote  sensing  can  sometimes  detect  agricultural 
diseases  and  pest  attacks  before  they  can  be  detected  from 
the  ground.  Bore  detailed  discussions  of  the  type  of 
instruments  used  and  examples  of  the  data  available  from 
remote  sensing  can  be  found  in  references  4,  6,  7,  13,  15, 
and  16. 

In  recent  years  land  has  emerged  as  the  one  resource 
"basic  to  both  the  good  life  and  the  goods  of  life"  [ref  d ]. 
There  has  been  emphasis  on  efficient  land  use  planning, 
particularly  in  metropolitan  areas.  Many  experts  are 
convinced  that  a well  conceived  program  of  remote  sensing  of 
metropolitan  areas  will  enhance  the  quality  of  urban  design. 
This  program  could  provide,  for  a given  area  of  land,  its 
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biological  pcoductivity,  soil  lioitations,  topology,  access 
to  transportation  and  utilities,  existing  vegetation  and  the 
use  of  adjacent  land  areas.  This  data  would  be  a very 
valuable  aid  in  efficient  urban  planning. 

NASA  has  been  involved  in  the  development  of 
instrumentation  to  perform  the  variety  of  sensing  reguired 
and  demonstration  of  the  data  available  to  potential  users. 
In  addition,  NASA  has  been  asked  by  potential  users  to 
estimate  total  program  cost  and  to  recommend  base  locations 
for  a special  fleet  of  aircraft  to  support  a program  of 
remote  sensing  of  urban  centers.  Assuming  an  initial  cost 
for  the  program  set  up  and  instrument  procurement  and  a 
fixed  cost  associated  with  each  sensing  mission,  NASA  has 
been  seeking  methods  of  analyzing  the  reduction  of  operating 
costs  through  efficient  vehicle  routing  and  effective  base 
locations. 

The  purpose  of  this  paper  is  to  develop  an  algorithm 
which  will,  given  a base  location  and  a set  of  urban  areas 
to  be  sensed,  provide  an  optimal  or  near  optimal  routing  for 
the  number  of  tours  reguired  to  visit  each  of  the  areas.  It 
is  assumed  that  the  operating  cost  will  be  primarily  a 
factor  of  distance  traveled  plus  a fixed  charge  for  each 
tour  required  to  visit  the  set  of  urban  areas.  Using  the 
fewest  flights  (subtours)  to  visit  each  area  and  optimally 
routing  the  aircraft  to  insure  the  minimum  distance  traveled 
will  yield  the  minimum  operational  cost. 

This  problem  was  formulated  as  a multiple  traveling 
salesman  problem  with  the  addditional  constraint  of  aircraft 
range.  The  aircraft  or  sensing  platform  may  only  be  capable 
of  visiting  a subset  of  the  set  of  urban  centers  before  it 
must  return  to  its  base  (completing  a subtour)  to  refuel. 
Thus  the  problem  is  to  minimize  the  sum  of  the  distances  of 
the  subtours.  The  algorithm  presented  by  Raymond  [ref  14] 
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was  modified  to  consider  this  range  constraint.  The 
modified  algorithm  was  programmed  in  FORTRAN  and  will  be 
discussed  along  with  operational  results  and  recommendations 
for  additional  uses. 
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II.  PHOBLEM  PORMOLATION 


The  traveling  salesman  problem  (TSP)  is  one  of  the 
classic  challenges  in  operations  research  (OR) . The  problem 
is  simply  stated:  a salesman  starting  from  his  home  wishes 
to  visit  a set  of  cities  once  and  only  once  and  return  home. 
In  what  order  should  the  cities  be  visited  to  minimize  the 
total  distance  traveled?  Mathematical  formulation  of  the 
problem  and  a discussion  of  the  constraints  can  be  found  in 
Wagner  [ref  17]  Chapter  13. 

OR  literature  has  included  applications  of  this  model  to 
a variety  of  practical  problems.  Industrial  organizations 
have  used  the  model  to^  improve  the  routing  of  delivery 
vehicles,  an  important  consideration  in  light  of  increased 
oil  prices  and  rising  driver  salaries.  The  "Chinese  postman 
problem"  is  a similiar  problem  applied  to  postal  service 
introduced  by  the  Chinese  mathematician  Kwan  [ref  10]. 
Municipal  governments  have  found  the  model  valuable  for 
routing  garbage  collection,  street  cleaning,  meter  reading 
and  snow  removal  vehicles.  Other  applications  include  dairy 
and  newspaper  delivery.  Further  discussion  of  applications 
can  be  found  in  Bradley's  "Survey  of  Deterministic  Networks" 
[ref  3]. 

The  key  difference  in  the  classic  statement  of  the  TSP 
and  the  remote  sensing  problem  is  the  addition  of  a range 
constraint.  The  salesman  of  the  TSP  may  refuel  during  his 
trip  without  altering  the  optimal  route.  This  is  not 
feasible  for  the  aircraft  performing  remote  sensing.  There 
exists  a point  during  the  trip  when  it  must  return  to  its 
nase  to  refuel  and  then  continue  the  mission.  Therefore, 
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the  problen  statement  becomes:  an  aircraft  starting  from  its 
base  must  visit  a set  of  urban  areas  once  and  only  once  and 
return  to  -its  base.  The  aircraft  is  subject  to  a range 
limitation  for  any  given  flight;  that  is,  it  may  only  be 
capable  of  visiting  a subset  of  the  set  of  cities  b&fore  it 
must  return  to  its  ba-se  to  refuel.  Each  time  the  aircraft 
returns  to  its  base  a subtour  is  completed  and  additional 
flights  (subtours)  are  created  as  necessary  to  accomplish  a 
visit  to  each  city.  The  objective  is  to  minimize  total 
distance  traveled  as  the  sum  of  the  length  of  each  subtour. 

A subtour,  t,  can  be  represented  as  a set  of  ordered 
city  pairs,  e.g. 

t = [ (i  ,i  ) , (i  ,i  ) . . . (i  ,i  ) , (i  ,i  ) ] 

1 2 2 3 n-1  n n 1 

where  i is  always  the  base.  The  mathematical  formulation 
1 

of  the  problem  is  to  find  the  set  of  subtours  to: 


minimize  \ 1 2.1 

k=1 


subject  to 


2-2 


1 < r 

k 


k~1 ,2,. . . ,s 


2.3 


1 < s < nc  2.4 

each  city  in  exactly  one  suhtour 
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where 


s s number  of  subtours 


1 = length  of  subtour  k 


n = number  of  cities  in  subtour  k 
k 


d = distance  from  city  i to  city  j 
ij 

r = range  of  aircraft 

nc  = number  of  cities  to  be  visited 

Note  that  if  s=1,  eg.  2.2  replaces  eg.  2.1  and  the  problem 
is  a standard  TSP.  The  upper  bound  of  eg.  2.4  represents 
the  worse  case  situation  that  each  city  oust  be  visited  by  a 
seperate  subtour.  Eg.  2.3  represents  the  range  constraint. 

The  distance  matrix  represents  the  euclidean 

distance  from  city  i to  city  j.  During  actual  sensing 

missions  the  aircraft  must  fly  several  passes  over  a city  to 

completely  image  the  area.  This  distance  (loiter  distance) 

depends  on  three  factors:  the  size  of  the  area,  the  percent 

side  overlap  used  for  film  images  and  the  turning  radius  of 

the  aircraft.  Bach  city  will  have  a unigue  loiter  distance 

(Id  ) that  must  be  considered  when  the  subtour  length  is 
i 

calculated.  Hhen  the  problem  is  formulated  the  value  Id  is 

i 

added  to  each  element  in  column  i.  Thus  the  new  distance 

matrix  is  D*  = D + Id  for  all  i.  The  new  matrix  is 
ij  ij  i 

used  to  compute  subtour  lengths. 


This  problem  appears  similiar  to  two  existing  expansions 


of  the  TSP;  the  aultisalesman  problem  (HTSP)  and  the  vehicle 
dispatching  problem.  In  the  MTSP,  m salesmen  starting  from 
a base  city  must  visit  n cities  and  return.  Each  city  is 
visited  by  one  and  only  one  of  the  m salemen.  The  objective 
is  to  route  each  of  the  m salesmen  such  that  the  total 
distance  traveled  is  minimized.  Bellmore  and  Hong  [ref  2] 
provide  a detailed  discussion  and  formulation  of  this 
problem.  The  objective  of  the  vehicle  dispatching  problem  is 
to  obtain  delivery  routes  from  a depot  to  a set  of  demand 
points  for  a fleet  of  vehicles  such  that  the  total  distance 
traveled  by  the  entire  fleet  is  minimized  as  with  the  m 
traveling  salesmen.  The  problem  is  enhanced  by  finite 
vehicle  capacities,  route  time  constraints  and  varying 
quantity  demanded  at  each  demand  point.  A complete 
description  cf  this  type  application  can  be  found  in  Golden 
[ref  12]. 

There  are  important  differences  between  the  remote 
sensing  problem  and  the  applications  mentioned  above.  The 
HTSP  has  no  upper  bound  on  the  number  of  cities  any  one  of 
the  m salesmen  can  visit.  The  aircraft  range  constraint 
requires  consideration  of  a upper  bound  on  the  number  of 
cities  that  can  be  visited  during  any  given  flight.  The 
vehicle  dispatching  problem  is  concerned  with  routing  each 
vehicle  in  the  fleet.  Although  potential  users  will 
presumably  have  a "fleet”  of  aircraft  to  accomplish  its 
sensing  missions,  not  every  aircraft  will  be  scheduled  to 
sense  a subset  of  the  areas  to  be  visited.  For  example, 
given  a fleet  of  three  aircraft  and  the  determination  that 
two  subtours  are  required  to  visit  the  set  of  cities,  only 
two  aircraft  (or  one  aircraft  making  two  flights)  will  be 
scheduled.  This  decision  assumes  the  fixed  cost  charged 
each  time  an  aircraft  is  launched  is  sufficiently  high  to 
make  it  more  desireable  to  create  the  fewest  number  of 
flights. 
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The  reoote  sensing  problem,  therefore,  is  not  a simple 
I restatement  of  an  existing  TSP.  The  addition  of  a range 
< constraint  requires  special  considerations.  The  solution 
1 technique  must  carefully  create  only  the  number  of  subtours 

required  and  route  each  subtour  so  that  the  total  distance 
traveled  is  minimized. 


III.  SOLOTION  ALGORITHMS 


j. 


One  o£  the  most  intriguing  aspects  of  the  ISP  is,  that 
although  siaply  stated,  it  renains  unsolved  in  a closed 
forn.  The  conputational  difficulty  with  known  solution 
tecnnigues  is  the  exponential  increasing  solution  tiae  as 
the  size  of  the  problem  increases.  There  are  two  types  of 
solution  algorithms:  heuristic  and  exact.  The  latter 
yields  the  optimal  tour.  The  difficulty  with  this  type 
algorithm  is  the  large  computer  storage  and  tiae 
requirements  which  are  related  to  the  large  number  of 
possible  solutions.  For  example,  given  a symmetric  TSP 
there  are  (n-1)!/2  possible  solutions  (n=nuaber  of  nodes).  A 
relatively  small  problem  of  ten  nodes  yields  181,440 
solutions.  Fortunately,  many  applied  or  real  world  problems 
do  not  demand  the  optimal  solution.  Powerful  heuristic 
algorithms  have  been  developed  which  will  yield  near  optimal 
solutions  and  provide  a substantial  savings  in  computer 
storage,  tiae  and  money.  An  excellent  survey  of  existing 
solutions  is  given  by  Bellmore  and  Nemhauser  [ref  2]. 

The  Bellmore  and  Nemhauser  article  classifies  the 
algorithms  by  solution  generation  method.  There  are  three 
fundamentally  different  ways  of  generating  solutions: 
tour-to-tour  improvement,  tour  building,  and  subtour 
elimination.  The  remote  sensing  problem  requires  tour 
length  feasibility  be  maintained  as  the  tour  is  formed. 
That  is,  if  visiting  the  next  city  (node)  will  result  in 
violating  the  range  constraint,  a new  subtour  is  started. 
Consequently,  a tour  building  solution  was  considered  the 
best  technique  to  adapt  to  this  problem. 
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Both  exact  and  heuristic  algocithns  were  studied  for 
applicability.  A literature  survey  provided  a number  of 
both  types  of  algorithms.  An  exact  solution  considered  was 
presented  by  Christofides  and  Eilon  [ref  5].  Their  paper 
discusses  three  solution  methods  to  the  vehicle-dispatching 
problem.  The  branch  and  bound  approach , a modification  of 
the  algorithm  published  by  Little  et.  al.  [ref  1 1 ]>  would 
reach  an  exact  solution  to  this  problem  provided  the  number 
of  subtours  required  was  known  a priori.  If  s (the  number 
of  subtours)  was  known,  the  base  location  would  be  replaced 
by  s artificial  bases  as  explained  in  ref  5.  The  algorithm 
would  progress  using  the  branch  and  bound  method  to  insert 
nodes  into  one  of  the  tours.  Before  branching  to  a new  node, 
it  would  be  necessary  to  check  that  the  total  distance 
accumulated  on  each  subtour  did  not  exceed  the  range.  The 
key  to  this  algorithm  is  knowing  the  number  of  subtours  a 
priori. 

As  noted  in  the  problem  formulation  section,  there 
exists  an  upper  bound  on  the  number  of  subtours.  Therefore, 
it  is  feasible  to  select  an  s known,  a priori  to  be 
sufficiently  large,  and  apply  this  branch  and  bound 
algorithm.  The  solution  would  be  recorded,  the  number  of 
subtours  decreased  by  one  and  a new  soluton  computed.  This 
procedure  would  continue  until  there  are  insufficient 
subtours  to  complete  a visit  to  each  city  and  the  last 
recorded  solution  would  be  the  optimal.  The  success  of  this 
procedure  would  be  very  dependent  upon  the  initial  selection 
of  s.  A poor  selection  could  involve  many  iterations  and 
would  require  considerable  user  interaction  and  computer 
time.  This  would  be  particularly  undesireable  for  problems 
that  require  a solution  in  a reasonably  quick  time. 
(Problems  of  this  type  are  discussed  in  Section  VI) . 
Therefore,  application  of  this  algorithm  was  not  considered 
profitable. 
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Liaiting  the  study  of  heuristic  algorithas  to  those 
using  the  tour  building  generation  technique,  the  algoritha 
presented  by  fiayaond  [ref  14]  was  selected  to  aodify  and 
apply  to  this  problea.  The  approach  of  Rayaonds's  algoritha 
is  siailiar  to  that  of  Karp  and  Xhoapson  [ref  9]  and  will 
obtain  optiaal  or  near  optiaal  solutions  to  the  classical 
TSP.  The  algoritha  successively  inserts  nodes  into  the  tour 
by  using  a series  of  decision  rules  and  provides  for  an 
exchange  procedure  to  deteraine  if  the  total  distance 
traveled  aay  be  reduced  by  rearranging  the  order  in  which 
the  nodes  are  visited.  k detailed  explanation  of  the 
algoritha  steps  is  available  in  ref  14. 

The  aodified  algoritha  starts  by  creating  a subtour  froa 

the  base  to  the  most  distant  city.  The  saae  criteria  as  the 

original  algoritha  is  used  to  select  the  next  node  for 

insertion.  That  is,  increment  calculations  (1  ) are  made 

Ic 

for  each  remaining  node  (k) . This  calculation  is  aade  by 

inserting  the  node  k between  each  pair  of  nodes  (p,q)  in  the 

subtour  and  computing  the  increaental  increase  in  subtour 

length,  I==d  -«-d  -d.  Each  node  will  have  a ainiaua 

k pk  kg  pq 

increment  calculation  (I*  ) and  a next  best  increment 

k 

calculation  (J*  ) . If  the  number  of  cities  in  the  tour  is 

k 

less  than  four  the  aaxiaua  I*  determines  the  node  k to  be 

k 

inserted  between  the  nodes  that  produced  I*^.  Otherwise, 

the  node  k for  which  J*  ->’1*  is  a aaxiaua  is  selected  and 

k k 

inserted  between  the  nodes  that  produced  Before 

completing  the  insertion,  the  algorithm  looks  ahead  to 

deteraine  if  the  range  constraint  will  be  violated.  If  it 
is  not  violated  the  node  is  inserted  and  the  three  link 
exchange  procedure  described  by  Raymond  is  applied.  This 
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1 


17 


i 


procedure  takes  each  node  in  turn  out  of  its  present 
subtour « rsevaluates  its  insertion  increaents  in  that 
subtour  and  inserts  the  node  between  the  nodes  that  now  give 
the  siniBUB  tour  length  increment.  If  the  constraint  is 
violated  the  node  is  not  inserted  euid  an  additional  subtour 
is  foraed. 

The  nearest  city  to  the  base  not  yet  in  a subtour  is 
used  to  fora  the  next  subtour.  The  criteria  for  inserting 
nodes  for  the  aultiple  tour  problem  is  changed.  A node  is 
inserted  in  the  subtour  that  will  create  the  least  increase 
in  total  distance  (sub  of  subtour  lengths) . A procedure  was 
developed  that  exchanges  nodes  aaong  subtours  to  deternine 

if  a savings  in  total  distance  traveled  is  possible  by 

• 

placing  a node  in  a different  subtour.  Before  any  node 
insertions  or  exchanges  aaong  subtours  is  coapleted,  the 
algoritha  looks  ahead  to  insure  tour  length  reasibility  is 
maintained.  The  algoritha  progresses  iteratively  until  each 
city  is  visited. 
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IV.  lai  ALGORITHM 


The  steps  of  the  algocitha: 

1.  With  the  base  as  one  node  select  the  node  farthest  from 
base  (assuaed  to  be  j*)  and  fora  a tour  with  lints 
and  Set  the  nuaber  of  subtours  s equal  to  1 and  the 

Quaber  of  lints  1 equal  to  2. 


2.  Define  I*  as  the  ainiaua  tour-length  increaent  that 
t 

results  froa  the  insertion  of  node  t in  each  existing  lint 

(P»*l)  • Define  nodes  p*  and  q*  as  the  nodes  connected  by  the 

lint  (p*#<3*)  that  produced  I*  . In  addition  define  J*  as 

t t 

the  increaent  for  the  node  t aost  nearly  equal  to  I*  . Thus 

t 

for  node  t,  I*  is  the  best  increaent,  J*  the  next  best 

t t 

increaent  and  (p*,q*)  the  lint  to  breat  if  inserting  node  t. 

Calculate  I*  and  J*  for  each  node  that  is  not  currently 
t t 

included  in  the  tour. 


3.  If  1 ^ 4 choose  the  node  a for  which  I*  is  a aaxiaua. 

1 t 

If  1 >4  choose  the  node  a for  which  J*  - I*  is  a 

1 t t 

aaxiaua. 


4.  Insert  the  node  a into  the  tour  between  p*  and  q*  by 
adding  lints  (p*,a)  and  (a,q*)  and  deleting  lint  (p*,q*)  . 


f 


Calculate  the  new  tour  distance  d^.  If  d^  S range,  coaplete 
the  insertion  and  increase  the  nuaber  of  links  in  the  tour 
by  one.  If  d^  > range,  do  not  insert  node  ■ and  go  to  step 
7 to  create  a new  subtour. 


5.  Take  each  node  in  turn  out  of  the  present  tour, 
reevaluate  its  insertion  increeents  and  insert  the  node 
between  the  pair  of  nodes  that  now  gives  the  ainiaua  tour 
length  increaent.  (In  aost  cases  the  node  is  reinserted  in 
the  link  created  when  the  node  was  taken  out  of  the  tour, 
producing  no  change.) 

6.  If  all  nodes  have  been  added  to  the  tour  stop,  otherwise 
go  to  step  2. 

7.  Select  the  node  a*,  not  already  in  a subtour,  closest  to 

the  base  and  fora  a tour  with  links  (B,a*)  and  (ae,B).  Set 

the  nuaber  of  subtours  s equal  to  s-»>1  and  the  nuaber  of 

links  1 equal  to  2. 
s 


6.  Calculate  I*^  as  in  step  2 for  each  node  k reaaining  to 
be  added  to  a subtour.  Consider  the  subtours  in  the  order 

in  which  they  were  created.  For  each  subtour  select  the 

node  that  produced  the  ainiaua  I*  . Before  inserting  the 

k 

node  two  checks  are  aade.  First  insure  that  the  range 

constraint  will  not  be  violated.  Second  insure  that  the 
total  distance  traveled  (sua  of  subtour  lengths)  will  not  be 
less  if  the  selected  node  is  inserted  in  a different 
subtour.  Only  those  nodes  which  pass  these  checks  are 
inserted.  Thus  at  any  given  iteration  a node  need  not  be 
inserted  in  each  subtour.  A given  node  aay  have  a saaller 
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I*  in  one  subtoar,  but  if  inserted,  would  force  range 
k 

constraint  violation.  In  this  situation  place  the  node  in 

the  next  best  subtour.  If  nodes  aay  not  be  added  to  any 
existing  subtour  without  violating  the  range  constraint  go 
to  step  7. 

9.  Bepeat  step  5 for  each  subtour  to  which  a node  was  added. 

10.  For  each  subtour  nake  the  increaent  calculations  using 
the  nodes  froa  the  other  subtours.  Select  the  ainiaua  I* 

k 

and  check  if  the  total  distance  traveled  nay  be  reduced  by 

adding  the  selected  node  to  the  given  subtour  and  deleting 
it  froa  its  original  subtour.  If  the  distance  is  reduced 
aake  the  change  and  repeat  step  5. 

11.  If  all  nodes  are  contained  in  a subtour  stop,  otherwise 
go  to  step  8. 
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TEST  PROBLEMS 


The  coaputer  progcaa  vas  used  to  provide  routes  froo 
bases  at  Moffett  Field,  California,  Kansas  City  and  Wallops 
Station,  Maryland  to  fifty  selected  urban  centers.  Figure  1 
is  the  grid  system  used  to  calculate  the  distances  between 
cities.  The  cities  were  placed  in  one  of  three  sections  by 
the  following  criteria:  the  straight  lines  from  Moffett 
Field  to  Kansas  City  and  Wallops  station  to  Kansas  city  were 
bisected  by  a perpendicular,  as  shown  in  Figure  1,  creating 
three  sections.  The  cities  would  be  imaged  from  the  base  in 
their  section.  Table  1 provides  the  coordinates  and  loiter 
distance  of  each  city.  The  range  of  the  aircraft  was 
assumed  to  be  2900  statute  miles.  Tables  2-1,  2-2  and  2-3 
are  the  routing  recommendations  of  the  computer  program. 

Tables  3-1,  3-2  and  3-3  are  the  results  of  a feasible 
solution  to  the  same  problem  computed  at  NASA's  Ames 
Research  Center.  The  constraints  observed  for  determining 
this  solution  were:  the  aircraft  would  have  a range  of  2900 
statute  miles,  no  more  than  three  urban  centers  would  be 
imaged  during  any  subtour  and  flights  from  Moffett  Field  and 
Wallops  Station  would  be  held  to  a reasonable  minimum. 

The  NASA  solution  required  29  subtours  and  a total 
distance  of  52,660  statute  miles.  The  algorithm  solution 
constructed  15  subtours  with  a total  distance  of  36,505 
statute  miles  traveled.  A savings  of  16,155  statute  miles 
results  from  the  algorithm  solution.  This  large  savings  is 
due,  in  part,  to  the  differing  criteria  used  to  formulate 
the  solutions.  However,  it  does  show  the  usefulness  of  the 
algorithm  and  the  speed  at  which  the  program  can  compute  a 
final  solution. 
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TABLE  1 


< CITY  COOEDINATES 

AND 

LOITEE 

DISTANCE 

City 

X 

Y 

Loiter 

Distance 

Albany 

2230 

1390 

190 

Allentown 

2190 

1230 

190 

Atlanta 

1950 

700 

290 

Baltimore 

2190 

1150 

340 

Birmingham 

1820 

650 

340 

Boston 

2350 

1470 

280 

Buffalo 

2030 

1310 

240 

Chicago 

1660 

1120 

630 

Cincinnati 

1830 

1015 

350 

Cleveland 

1930 

1160 

410 

1 Columbus 

1900 

1080 

230 

: Dallas 

1350 

500 

740 

■*  Dayton 

1820 

1080 

250 

Denver 

I 

910 

900 

450 

Detroit 

1850 

1230 

350 

Greensboro 

2150 

910 

290 

Hartford 

2300 

1350 

150 

Houston 

1470 

350 

670 

* Indianapolis 

1730 

1030 

400 

• Kansas  City 

1380 

880 

380 

i Los  Angeles 

220 

560 

670 

I Louisville 

« 

1780 

910 

190 

1 

j Hemphis 

1560 

710 

220 

1 Miami 

2300 

3 70 

370 

1 Milwaukee 

1 

1620 

1210 

270 

! 1 

Loiter  distances  provide  by 

Ames 

Research  Center. 
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TABLE 

1 (con* 

d) 

CITY  COOHDINATES  AND 

LOITEB 

DISTANCE 

City 

X 

Y 

Loiter 

Distance 

Hianea polls 

1400 

1260 

380 

Nashville 

1780 

800 

270 

New  Orleans 

1740 

410 

290 

New  York 

2300 

1300 

730 

Norfolk 

2270 

2020 

230 

Oklahoma  City 

1290 

660 

330 

Omaha 

1300 

1010 

230 

Philadelphia 

2230 

1200 

480 

Phoenix 

560 

4 70 

250 

Pi  ttsburgh 

2030 

1150 

430 

i Portland 

150 

1250 

450 

Providence 

2360 

1370 

160 

Bochester 

2090 

1330 

150 

. Sacramento 

150 

850 

410 

I 

St  Louis 

1590 

910 

500 

St  Petersburg 

2140 

410 

250 

Salt  Lake  City 

5 80 

910 

200 

San  Antonio 

1270 

270 

290 

* San  Bernardino 

310 

560 

190 

1 San  Diego 

320 

480 

510 

i San  Francisco 

100 

7 80 

530 

i Seattle 

1 

190 

1420 

500 

i Syracuse 

2160 

1350 

310 

1 Toledo 

1830 

1150 

230 

1 Washington  DC 

2200 

1100 

350 

Moffett  Field 

120 

760 

Kansas  City 

1380 

880 

' Wallops  Station 

2270 

1080 

1; 


/' 


j. 


I 
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TABLE  2-1 

FLIGHT  Missions  POH  ORflAN  CENTER  COVERAGE 
BASE:  MOFFETT  FIELD 

Subtour  Cities  Imaged 

1 Sacramento  - Seattle  - Portland 

2 Los  Angeles  - San  Diego  - San  Bernadino 

3 San  Francisco  - Salt  Lake  City  - Phoenix 


i 

i Total  Distance  = 7772  miles 

I Time  to  Compute  Solution=  1.27  seconds 


Length 

2892 

2184 

2696 


! 

I 

j 

I 

•j 

\ 

i 

i 


) 


• J 
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TABLi:  2-2 


FLIGHT  HISSIONS  FOR  URBAN  CENTER  COVERAGE 
BASE:  KANSAS  CITY 


Sabtour 

1 

2 

3 

4 

5 

6 


Cities  Imaged 

Houston  - San  Antonio  - Oklahoma  City 

Dallas  - Denver 

Omaha 

Indianapolis  - Louisville  - St.  Louis  - Kansas  City 
Nashville  - Birmingham  - New  Orleans  - Memphis 
Minneapolis  - Milwaukee  - Chicago 


Total  Distance  = 14092 


Length 

2878 

2853 

581 

2520 

2746 

2514 


Time  to  Compute  Solution  = 6.97  seconds 


TABLE  2-3 


PLIGHT  HISSIONS  FOR  OBBAH  CENTEB  COVERAGE 
BASE:  WALLOPS  STATTION 


Subtour 

Cities  Imaged 

Length 

1 

St  Petecsberg  - Hiani 

- Norfolk 

2642 

2 

Columbus  - Cincinnati 

- Dayton  - Detroit 

- Buffalo 

2816 

3 

Hartford  - Providence 

- Boston  - Albany  - 

Philadelphia 

2401 

4 

Mew  YorJc  - Syracuse  - 

Rochester  - Allentown  - Baltimore 

2607 

5 

Washington,  DC  - Pittsburgh  - Cleveland  - 

Toledo 

2451 

6 

Greensboro  - Atlanta 

1724 

Total  Distance  = 14641 

Time  to  Compute  Solution  = 17.02  seconds 
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TABLE  3-1 

FLIGHT  HIS SIGHS  FOB  OBBAH  CENTEB  C07EBAGE 


BASE:  aOFFETT  FIELD 

Sobtoar 

Cities  Imaged 

LengtJ 

1 

Salt  Lake  City 

1950 

2 

Phoenix 

2340 

3 

Denver  - Oiaha 

1990 

4 

Houston 

2020 

5 

San  Antonio 

1990 

6 

Dallas 

1990 

7 

Oklahoma  City  - Kansas  City 

1410 

8 

Mew  Orleans  - Hemphis 

1960 

9 

St  Louis  - Minneapolis 

1850 

10 

Chicago  - Milwaukee 

1800 

11 

Miami 

2780 

12 

St  Petersburg 

2250 

13 

Birmingham  - Atlanta 

2070 

14 

Nashville  - Louisville  - Indianapolis 

1950 

15 

Cincinnati  - Dayton  - Columbus 

2100 

16 

Detroit  - Toledo 

1880 

17 

Cleveland  - Pittsburgh 

2380 

18 

Buffalo  - Bochester 

2190 

19 

Syracuse  - Albany 

2710 

Total  Distance  = 39520  niles 


TABLE  3-2 


I 

FLIGHT  HISSIOHS  FOR  UBBAM  CENTER  COVERAGE 
BASE:  KANSAS  CITY 


sabtour 

Cities  laaged 

Length 

1 

Seattle 

1950 

2 

Portland 

1450 

3 

San  Francisco  - Sacraaento 

1140 

4 

Los  Angeles  - San  Bernadino 

1510 

5 

San  Diego 

1400 

Total  Distance  » 7450  ailes 


TABLE  3-3 


PLIGHT  Missions  FOR  OHBAH  CENTER  COVERAGE 
BASE:  WALLOPS  STATION 


Subtour 

Cities  laaged 

Length 

1 

Boston  - Providence  - Hartford 

1430 

2 

New  York 

1120 

3 

Philadelphia  - Allentown 

1080 

4 

Washington,  DC  - Baltiaore 

970 

5 

Norfolk  - Greensboro 

1090 

Total  Distance  = 5690  alles 


VI . ADDITIONAL  APPLICATIONS 


The  purpose  of  this  paper  was  to  develop  an  algoritha 
that  will  provide  optiiaal  or  near  optimal  routing  for  a 
range  constrained  vehicle  visiting  a set  of  nodes.  The 
algorithm  may  be  used  to  solve  similiar  problems.  Several 
considerations  for  use  of  this  algorithm  are  discussed  in 
this  section.  They  are:  under-flying  satellite  coverage, 
isolated  disaster  coverage  and  point  reconaissance . 

Satellites  are  often  used  to  image  large  land  masses  and 
bodies  of  water.  These  orbiting  platforms  perform  the 
sensing  mission  during  periodic  passes  over  the  critical 
area.  The  product  of  a mission  may  sometimes  be  impaired  by 
cloud  cover  or  may  reveal  the  need  for  finer  resolution 
imagery.  The  areas  requiring  additional  coverage  are 
assumed  dispersed  and  relatively  small.  These  areas  can  be 
considered  as  the  cities  of  the  original  problem  and  the 
algorithm  used  to  determine,  from  a given  base,  the  number 
of  flights  required,  the  optimal  or  near  optimal  routing  and 
the  total  distance  traveled. 


I Disasters  such  as  tornados  or  flooding  rivers  may 
I require  government  agencies  to  estimate  damage  and  make 
; recommendations  for  appropriating  the  relief  effort.  Remote 
I sensing  can  provide  valuable  data  to  assist  in  the  decision 
! process.  consider  a flooding  river;  it  is  likely  several 
I areas  along  the  flood  path  need  immediate  attention.  The 

goal  is  to  provide  the  data  for  evaluation  as  rapidly  as 
possible.  The  algorithm  can  provide  routing  from  a base  to 
the  areas  of  concern  such  that  the  sensing  mission  is 
completed  in  the  shortest  possible  time.  It  is  assumed  that 
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. I 

« 

least  distance  routing  iaplies  the  least  tiae  will  be 
I required  to  accoeplish  the  aission. 

Several  other  applications  are  sioiliar  to  the  disaster 
coverage  and  involve  routing  reconnaissance  aircraft.  For 
exaaple,  when  fighting  large  forest  fires,  "hot  spots'*  | 

develop  that  need  periodic  aonitoring.  Also,  in  a large  j 

theater  of  operations,  the  ailitary  coanander  often  requires  i 

tiaely  reconneissance  of  key  objective  areas.  In  both  these  | 

situations  the  goal  is  to  provide  the  inforaation  to  the 
user  as  rapidly  as  possible.  Therefore,  the  algoritha  can 
be  used  to  route  the  aircraft  to  the  critical  areas  and  | 

provide  the  inforaation  as  quickly  as  possible. 
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VII.  COMCLDSIONS 


The  algoritho  aad  computer  program  presented  can 
successfully  construct  an  optimal  oc  near  optimal  solution 
to  the  range  constrained  vehicle  routing  problem  formulated 
in  this  paper. 

Several  aspects  of  remote  sensing,  considered  beyond  the 
scope  of  this  paper,  provide  a basis  for  further  study  in 
this  area.  The  angle  of  the  sun  can  have  an  appreciable 

effect  on  the  resolution  of  the  imagery.  Sun-angle  is  a 
function  of  the  time  of  day  and  could  be  considered  in  the 
computer  program  by  adding  a time  dimension.  Using  aircraft 
speed,  it  would  be  possible  to  calculate  a time  period 
during  which  the  aircraft  would  image  a given  city.  This 
would  provide  the  user  the  ability  to  insure  that  the 
desired  sun-angle  occurs  during  the  imaging. 

This  paper  assumes  each  aircraft  has  sufficient  sensor 

capacity  to  image  any  number  of  cities  in  a subtour.  It  is 

likely  there  may  be  a finite  capacity  for  certain  sensors. 

For  example,  assume  only  a fixed  amount  of  film  can  be 

I placed  on  the  aircraft.  The  film  reguied  to  image  all  the 

cities  in  a given  subtour  may  exceed  this  amount.  To 

i consider  this  problem,  the  aircraft  would  start  with  a 

I finite  amount  of  film  which  would  be  reduced  by  the  amount 

required  to  image  a given  city  as  that  city  is  visited. 

j Thus  sensor  capacity  could  be  added  as  a constraint  and 

i considered  in  the  same  manner  as  the  range  constraint. 

i 

I 

i The  problem  of  imaging  areas  as  widely  dispersed  as  the 

I fifty  urban  centers  listed  in  Table  1 requires  multiple 

I 

I 
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bases.  The  example  pcoblem  assigned  cities  to  bases  by 
section.  This  approach  is  reasonable,  however,  it  would  be 
preferable  to  permit  the  computer  program  to  make  this 
assignment.  One  consideration  is  to  construct  an  artificial 
base  from  which  all  flights  originate  and  sufficient  dummy 
bases  at  the  sane  location  as  the  actual  bases  to  permit 
multiple  subtours  through  that  base.  The  distance  matrix 
must  be  constructed  such  that  each  subtour  starts  from  the 
artificial  base  to  one  of  the  actual  bases  and  is  routed 
back  through  the  actual  base  to  the  artificial  base. 
Formulation  of  a fifty  city  problem  with  artificial  and 
dummy  bases  could  require  considerable  computer  storage. 
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APPENDIX  A 


COMPUTER  PROGRAM 


variable  definitions 
xm»x  coordinate  of  the  i-th  city 

Y(I)»Y  COORDINATE  OF  THE  I-TH  CITY 
Zm*LOITER  DISTANCE  OF  THE  I-TH  CITY 
0(1* J)»DISTANCE  MATRIX 
RANGE»AIR.CRAFT  RANGE 
NTOURS=NUMBER  OF  SUBTOURS  CREATED 
LINKS!  n»NJMBER  OF  LINKS  IN  SUBTOUR  I 
INNODE (I, J)*VECTOR  OF  CITIES  IN  SUBTOUR  I 
NONODE(I)»VECrOR  OF  CITIES  TO  BE  VISITED 
LENGTH! I )»LENGTH  OF  SUBTOUR  I 
TOTDIS»TOTAL  DISTANCE  TRAVELED 

LH!I*J»*LEFT  HAND  ENDPOINT  OF  SUBTOUR  I*  LINK  J 

RH!I , JI=RIGHT  HAND  ENDPOINT  OF  SUBTOUR  I*  LINK  J 

NC=NUMB£R  OF  CITIES  IN  THE  PROBLEM 

S!I,J)*IN  SUBTOUR  I THE  SUCCESSOR  OF  CITY  J 

PRED!I,J)=IN  SUBTOUR  I THE  PREDECESSOR  OF  CITY  J 

> 

> 

> 

> STORAGE  MATRICES  FOR  ALGORITHM 

> calculations 

> 

> 


STAR!! 
NNSTR! t , J) 
LLSTR!I,J) 
SECSTR!! ,J) 
SNNSTR!I,J) 
SLLSTRd  , JI 
THETA! I, J) 


MAIN 


COMMON  0(25*25) *S!25* 25). I NNODE! 25*25 ) *N0N0DE!25) * 
1RH!25*25) *LH! 25* 25 ,L INKS! 25 ) » THATA!25 *25 ) »NTOURS.KK* 
1IEND*NC,TQT0IS,LENGTH!25) .RANGE* I STOP,  PREO ( 25 , 25 ) 
DIMENSION  X!25) ,Y!25) ,STAR!25*25) ,NNSTR!25*25), 
ILLSTR!25*25) , SECSTR(25*25)*SNNSTR(25,25) , SLLSTR! 25*25) 
1,CRIT!25) *0LDRH(25*25) *OLDLH( 25*25 )*Z! 25) 

INTEGER  S,RH*  OLDRH,OLOLH*  EXNOOE  * PREO 
REAL  LENGTH 
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c 


c 

c 


c 


c 

c 

c 


c 

c 

c 

c 


INITIALIZE  ALGORITHM  VARIABLES 


CALL  INIT 
IPRINT=0 
IFIISTOP  .NE. 


D)  GO  TO  100 


THE  CRITERIA  TO  FORM  THE  FIRST  SUBTOUR  IS  TO  SELECT  THE 
CITY  furthest  from  THE  BASE  AS  THE  FIRST  CITY  TO  VISIT 

AMAX=-1.0 
00  17  J=1,NC 

IF(D{1,J)  .3T.  AMAX)  GO  TO  18 
GO  TO  17 

18  AMAXsOdtJ) 

NOOE-J 
17  CONTINUE 

FORM  THE  FIRST  SJBTOUR 

S(  ItU^NOOE 
S(1,N00E)=*1 
LINKS! 1)*2 

ESTABLISH  ENDPOINTS 

53  CALL  ENOPTS 

IF  USER  DESIRES  TO  PRINT 
ITERATION  SET  IPRINT=3 


SUCCESSOR  FUNCTION  AT  EACH 


IFilPRINT  .NE.  3)  GO  TO  5^ 

00  109  I=»1,NT0URS 

WRITE(6»108)  (S(  If J),J-1»NC) 

108  F0RMAT(25I3) 

109  CONTINUE 

CALCULATE  TOUR  LENGTHS  AND  TOTAL  DISTANCE 
54  CALL  TORDIS 

ESTABLISH  NONOOE  AND  INNOOE  VECTORS 
CALL  TOVIST 

WHEN  IEND=0  ALL  CITIES  HAVE  BEEN  VISITED 
IFdENO  .EQ.  0)  GO  TO  50 
SELECT  criteria  TO  DETERMINE  NEXT  CITY  TO  VISIT 


999  STOP 

IFINTOJRS  .GT. 
IF(NTOURS  .EQ. 
CALL  ENOPTS 
CALL  THREEX 
GO  TO  53 


1)  CALL  GRATOR 
1)  CALL  SNGLTR 


100  WRITE(b,105)  ISTOP 

105  FORMAT! lOXf 'CITY* , I5t*  IS  OT  OF  RANGE*) 

GO  TO  999 

50  DO  110  I*lfNTOJRS 

WRITE(6fILl)  (Sdf  J)  fJ-l«NC) 

111  F0RMAT(25I3) 

110  CONTINUE 

DO  112  I=1,NT0JRS 

WRITE(6.113)  IfLENGTHd) 

113  FORMAT! lOX, ‘LENGTH  OF  SUBTOUR 13 t • *‘,F15.4) 

112  CONTINUE 
WRITE!6,115)  TOTDIS 

115  FORMAT! lOX, ‘TOTAL  01  STANCE ‘ , F20 .4 ) 

END 
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SUBROUTINE  INIT 


1 J 


C 

c 

c 


c 

c 


c 

c 


SUBROUTINE  WILL  ^EAO  PROBLEM  SIZE»  CITY  LOCATIONS  AND 
LOITER  DISTANCES,  AIRCRAFT  RANGE,  INITIALIZE  SUCCESSOR 
AND  PREDECESSOR  FUNCTIONS  AND  COMPUTE  DISTANCE  MATRIX 


COMMON  D( 25,25) ,S (25,25 ), INNOOE (25,25) ,NONODE( 25) , 
IRH(25,25) ,LH( 25,25.LINKS(25) , THATA(25,25) , N TOURS. KK, 

II  END, NC.TOTOIS. LENGTH (25) .RANGE, I STOP,PRED( 25,25) 
DIMENSION  X(25) ,Y(25),STAR(25,25) .NNSTR ( 25 , 25 ) , 
1LLSTR(25,25) , S= CSTR ( 25 ,25 ) , SNNSTR ( 25 , 25 ) , SL LSTR( 25 , 25 ) 
l,CRIT(25),0LDRH(25,25) ,OLDLH(  25, 25), Z(  25) 

INTEGER  S,RH,OLDRH,OLOLH, EXNODE, PRED 
REAL  LENGTH 


IPRINT=«0 

ISTOP*0 

NTOURS*! 


READ  NUMBER  OF  CITIES  IN  PROBLEM  AND  AIRCRAFT  RANGE 


READ  (5,11)  NC, RANGE 
11  FORMAT!  I3,F20.<V) 
HALFR=RANGE/2. 


READ  CITY  LOCATIONS  AND  LOITER  DISTANCES,  BASE  MUST  BE 
READ  AS  CITY  1 


DO  20  I*l,NC 

RE  AO (5, 21)  X( I),Y(I  ),Z(I  ) 
21  FORMATOFIO.^) 

20  CONTINUE 


COMPUTE  DISTANCE  MATRIX 

THE  FACTOR  1.15  CONVERTS  NAUTICAL  MILES  TO  STATUTE  MILES 


DO  12  I»1,NC 
DO  13  J*l,NC 

0(1,  J)*(SaRT(( (X(I  )-X(J) )#*2+((Y(I)-Y(J)  )**2)) ))♦ 
l.l3+Z(J) 


CHECK  FOR  CITY  OUT  OF  RANGE 


IF(I.EQ.l)  30  TO  14 
GO  TO  13 

14  IF(D(I,J)  .GT.  HAlFR)  IST0P=J 
13  CONTINUE 
12  CONTINUE 


INITIALIZE  SUCCESSOR  AND  PREDECESSOR  FUNCTION 


DO  15  I»l,NC 
DO  16  J»1,NC 
S( I, J)»0 
PRED(I,J)»0 
16  CONTINUE 
15  CONTINUE 


IF  USER  DESIRES  TO  PRINT  DISTANCE  MATRIX  SET  IPRINT=1 

1)  GO  TO  30 


IFdPRINT  .NE. 

DO  150  lal.NC 

WRITE(6)106)  (0(1, J) ,J»1,NC) 
106  F0RMAT(25F4.0) 

150  CONTINUE 
30  RETURN 
END 


SUBROUTINE  ENOPTS 


38 


'TV 


C 

C 


SUBROUTINE  ESTABLISHES, 
ENDPOITS  OF  EACH  LINK 


FOR  EACH  SU6T0JR,  THE 


I 1 


COMMON  0( 25,25) ,5(25,25) , I NNODE ( 25 ,25 ) ,N0N0DE(25) , 
1RH(25 ,25) ,LH( 25, 25, LINKS (25) , TH AT A ( 25 , 2 5 ) ,NTOURS,KK, 
lI£NO,NC,TOTDIS,LENGTH(25) , RANGE, ISTOP,PRED( 25,25) 
DIMENSION  X(2  5) ,Y(25) , STAR (25,25)  ,NNSTR ( 25 , 25 ) , 
ILLSTR(25,25) .SECS TR( 25,25 ) ,SNNSTR( 25, 25) , SL LSTR ( 2 5, 25 ) 
1, CRI T( 25 ), OLD RH( 25,25) , OLDLH( 25 , 25 ) , Z ( 2 5 ) 

INTEGER  S, RH,OLDRH,OLDLH, EXNODE, PRED 
REAL  LENGTH 


IPRINT=0 


23 

22 


DO 

22  J 

=1,NT0URS 

LH(  J 

,1)  = 1 

RH(  J 

,1)=S( J, 1 

LIMI 

T=LINKS(U 

DO 

23  K 

=2, LIMIT 

IM=K 

-1 

IL=RH( J,IM) 

LH(  J 

,K)=RH(J, 

IR=LH( J,K) 

RH(  J 

,K)*S( J,  I 

CONTINUE 

CONTINUE 


EACH  LINK  AT 


C IF  USER  DESIRES  TO  PRINT  THE  ENDPOINTS  OF 
C EVERY  ITERATION  SET  IPRINT=A 

IFdPRINT  .NE.  <►)  GO  TO  25 
DO  110  I=1,NT0JRS 
LIMIT=LINKS( I ) 

WRITE (6,102)  (LH( I, J ),RH(I , J) ,J=1,LIMIT ) 

102  F0RMAT(5X, 'ENDPOINT* ,2110) 

110  CONTINUE 
25  RETURN 
END 

SUBROUTINE  TOROIS 

C SUBROUTINE  CALCULATES  LENGTH  OF  EACH  SUBTOUR  AND  TOTAL 

C DISTANCE  AS  SUM  OF  SUBTOUR  LENGTHS 

COMMON  D( 25,25) ,S( 25,25 ), INNODE ( 25, 25 ) , NONODE ( 25 ) , 
1RH(25 ,25)  ,LH(  25 , 25  ,L I NKS ( 25 ) , THATA ( 25 , 2 5 ) ,NTOURS, 
1IEND,NC,T0TDIS,LENGTH(25)  ,RANGE, IST0P,3RED( 25,25) 
DIMENSION  X(25),Y(25) , ST AR ( 25 , 25 ) ,NNSTR ( 25 , 25 ) , 

1LLSTR( 25,25) , S= CSTR ( 25 , 25 ) , SNNS TR ( 25 ,25 ) , S L LSTR ( 2 5 , 2 5 ) 
1, CR IT ( 25 ) ,OLDRH( 25,25 ),OLOLH( 25,25) ,2(25) 

INTEGER  S,RH,OLDRH,OLDLH,EXNODE,PRED 
REAL  LENGTH 

1PRINT=0 
T0TDIS=0.0 
DO  55  I=1,NT0URS 
NS=S( 1, 1) 

NP=1 

LENGTH( I)=0.0 

20  LENGTH( I )=0(NP,NS)+LENGTH( I ) 

NP*NS 
NS=S( I,NP) 

IF(NP  .GT,  1)  GO  TO  20 
TOTDIS=TOTOIS^-LENGTH(I ) 

55  CONTINUE 


C 

C 


IF  USER  DESIRES  TO  PRINT  LENGTH  OF  EACH  SUBTOUR  AND 
TOTAL  DISTANCE  AT  EVERY  ITERATION  SET  IPRINT=5 


19 


mamgm 


ooo 


( 


ll  IFCIPRINT  .NE.  5)  GO  TO  25 

ll  00  105  I=1,NTDJRS 

jl  WRITE{6,132)  I,LEMGTH(I) 

ll  102  FORMATdOX,  • LENGTH  OF  SU  BTOOR  • , 13  , • =',F15.4) 

f 105  CONTINUE 

il  WRITE(6,101)  TOTOIS 

101  FORMATdOX, 'TOTAL  01  ST  ANCE  • , F20 . 4 ) 

25  RETURN 

ENO 

SUBROUTINE  TO\/IST 

SUBROUTINE  ESTABLISHES  A VECTOR  INONOOE)  OF  THE  CITIES 
THAT  HAVE  NOT  BEEN  VISITEO  ANO  FOR  EaCH  SUBTOUR  A VECTOR 
(INNODE)  OF  CITIES  IN  THAT  SUBTOUR 


COMMON  0(25,25) ,S (25,25), INNODE( 25, 25),  NONODE (25) , 
1RH(25,25)  ,LH(  25 , 25 , L I NKS( 25 ) , THATA ( 25  ,2 5 ) ,NTOURS, 

II END,NC,TOTOiS,  LENGTH (25) , RANGE, ISTOP,PREO(  25,25) 
DIMENSION  X(25) ,V(25) , STAR ( 25 , 25 ) , NNSTR ( 25 , 25 ) , 

1 LLSTR ( 25, 25 ), SECS TR( 25,25) , SNNS TR ( 25 , 25 ) , SL LSTR ( 2 5 , 25 ) 
1,CRIT(25) ,OLORH( 25,25) , OL DL H ( 25 , 25 ) , Z ( 25  ) 

INTEGER  S,RH,OLORH,OLDLH,EXNOOE,PRED 
REAL  LENGTH 


I EN0=0 

DO  24  J=2,NC 
DO  25  I=1,NT0URS 

IF(S(I,J)  .NE.  0)  GO  TO  24 

IF(I  .LT.  NTOURS)  GO  TO  25 

IEND=1 

KK=KK+1 

NONODE(KX) =J 

25  CONTINUE 
24  CONTINUE 

DO  26  1=1, NTOURS 
MM  = 1 

DO  27  J=2,NC 

IF(S(I,J)  .EQ.  0)  GO  TO  27 
I.NNODE(I,MM)=J 
MM=MM+1 
27  CONTINUE 

26  CONTINUE 
RETURN 
END 

SUBROUTINE  SNGLTR 

C SUBROUTINE  PLACES  NEXT  CITY  IN  SUBTOUR 


COMMON  0( 25,25) ,S (25,25) , INNODE (25,25) ,N. NODE (25) , 
1RH(25,25) ,LH( 25, 25 ,L INKS{ 25 ) , THATA (25,2 5) , NTOURS, 
1IENO,NC,TOTOI S,LENGTH(25) .RANGE, I STOP , PRED ( 25 ,25 ) 
DIMENSION  X(25),Y(25),STAR( 25,25) ,NNSt' ( 25,25) , 

ILLSTR  (25,25)  , SEC  STR  ( 25 , 25  ) , SNNSTR(25,25  ) , SL  I,  STR  ( 2 5 , 25 ) 
1 ,CRIT( 2 5) ,OLDRH( 25,25) , OLDLH( 25 , 25 ) , Z ( 25  ) 

INTEGER  S,RH,OLDRH,OlDLH, EXNODE, PRED 
REAL  LENGTH 


CALL  TNCRMT 
AMAX=99999. 

AMIN=-1.0 

IPASS=0 

IF(LINKS( 1)  .3T.  4)  GO  TO  34 
LIMIT=LINKS(NTOJRS) 


28  DO  29  J=1,KK  ' 

DO  30  1=1, LIMIT  . i 

1 


31 


IF(THETA(I,J)  .LT.  AMAX)  G3  TO  31 
GO  TO  30 
NN*^^ONODE(  J) 

LL=I 

AMAX  = THETA( I , J) 

30  CONTINUE 

STARd,  J)  = AMAX 
NNSTR( l,J) =NN 
LLSTR(l.J) =LL 
AMAX=99999. 

29  CONTINUE 

DO  32  K=1,KK 

IF(STAR(1,<)  ,GT.  AMIN)  GO  TO  33 
GO  TO  32 

33  IN=NNSTR(1,<) 

INL=LLSTR( I ,K) 

AMIN=STAR(1,K) 

32  CONTINUE 

IFdPASS  .EQ.  1)  GO  TO  35 
CHECK=LENGTH(  D+AMIN 
IFICHECK  .GT.  RANGE)  CALL  TOURNU 
IFtCHECK  .GT.  RANGE)  GO  TO  53 
IL=LH(NTOURS, INL) 

IR=RH( NTOURS, INL) 

S(NTOURS,IL)=IN 
S (NTOURSt IN)=IR 
LINKSd  ) = LINKS11)+1 
GO  TO  53 


34  IPASS=L 
GO  TO  2 8 

35  LIMIT=LINKS(NT0URS) 

AMAX=99999. 

DO  36  1*1, KK 
DO  37  J=l, LIMIT 

IF(THETA(  I,J)  .GT. STAR!  1,1)  .AND.THETAd  ,J).LT.AMAX) 
1 GO  TO  38 
GO  TO  37 

38  AMAX  = THETA( I , J ) 

LNK=J 

NOD=NONODE(<K) 

37  CONTINUE 

SECSTRd,I)*AMAX 
SNNSTRd,!  )=NOD 
SLLSTRd,!  )*LNK 
AMAX*99999. 

36  CONTINUE 

AMIN=-1. 

DO  39  K=1,KK 

CRIT  (X)  = SECSTRd,X  )-STAR(  1,K) 

IF(CRIT(K)  .GT.  AMIN)  GO  TO  40 
GO  TO  39 

40  AMIN=CRIT(<) 

TEST=STARd,<) 

NEWLNK*LLSTRd  ,K) 

NEWNOD*NNSTRd,K) 

39  CONTINUE 
CHECK*LENGTH(1)+TEST 

IFCCHECK  .GT.  RANGE)  CALL  TOURNU 
IFICHECK  .GT.  RANGE)  GO  TO  53 
1L*LH( NTOURS, NEWLNK) 

IR*RH(NTOURS,NEWLNK) 

S(NTOURS,IL)=NEWNOD 
S(NTOURS,NE(<NOO)*IR 
LINKSd)*LINKS(l)+l 
53  RETURN 
END 

SUBROUTINE  GRATOR 

C SUBROUTINE  SELECTS  NEXT  CITY  OR  CITIES  TO  BE  VISITED 


hi 


C AND  places  them  IN  THE  APPROPRIATE  SUBTOJR 


rr 


ii 


.LT.  AMAX)  GO  TO  74 


COMMON  0(25*25) , S ( 25 .25 }. INNODE ( 25,25 ) tNONOOE (25) » 
1RH(25,25I ,LH(25,25.LINKS(25) ,THATA(25,25) .NTOURS, 
lIEN0,NCiTQTDIS,LENGTH(25) .RANGE. ISTOP,PRED( 25,25) 
DIMENSION  X(2  5) .Y(25) ,STAR(25.25)  ,NNSTR ( 25 .25 ) . 

I LLSTR( 25,25) .SECSTRl 25.25 ) . SNNSTRi 25 , 25 ) ,SLLSTR(25 
l.CRIT(25 ) .OLORH(25,25) . OLDLH ( 25 , 25 ) , Z( 25 ) 

INTEGER  S.RH.OLORH, OLDLH, EXNOOE.PRED 
REAL  LENGTH 

ISAVTR=NTOURS 
DO  71  N=1,ISAVTR 
NTOURS=N 
LIMIT=LINKS(N) 

CALL  TNCRMT 
AMAX =99999. 

DO  72  J=1,KK 
DO  73  1=1, LIMIT 
IF(THETA(  I , J ) 

GO  TO  73 

74  NN=NONODE(J) 

LL  = I 

AMAX=THETA(I, J) 

73  CONTINUE 

STAR  (N,  J)=AMAX 
NNSTR( N, J) =NN 
LLSTR(N.J)=LL 
AMAX=99999. 

72  CONTINUE 

71  CONTINUE 

NTOURS=ISAVTR 
IFLAG=0 

DO  75  1=1, NTOURS 
AMAX=99999. 

00  76  J=l,KK 
IF(STAR(I , J) 

GO  TO  76 

77  1N=NNSTR( I ,J) 

INL=LLSTR(I, J) 

AMAX  = STAR(  I. J) 

ICOMPR=J 

76  CONTINUE 

00  78  K=l, NTOURS 

IF(<  .EQ.  I)  GO  TO  78 
IF(AMAX  .LT,  STAR(K, ICOMPR) ) GO  TO  78 
CHEC K=STAR ( < , I COMPR ) +LENGTH( X ) 

IF(CHECK  .LT.  RANGE)  GO  TO  75 

78  CONTINUE 

CHECK=AMAX«-LENGTH(  I ) 

IF  (CHECK  .GT.  RANGE)  30  TO  79 
IFLAG=1 
IL=LH(  I ,INL) 

IR=RH(I,INL) 

S( I,IL)=IN 
S(  I,  IN)  = IR 
LINKS( I )=LINKS( I)+l 
GO  TO  75 

79  STAR( I.IC0MPR)=99999. 

75  CONTINUE 

IFdFLAG  .EQ.  3)  CALL  TOJRNU 
CALL  ENOPTS 
CALL  TOROIS 
CALL  TOVIST 
CALL  CRSEXC 
RETURN 
END 

SUBROUTINE  TN:RMT 


LT.  AMAX)  GO  TO  77 


):2 


.25) 


i 

J 


ooo  no  non 


SUBROUTINE 
OF  PLACING 
sue  TOUR 


CO'I’JTES  THE  INCRE'iENTAL  DISTANCE  INCREASE 
ANf  GIVEN  NODE  INTO  ANY  GIVEN  LINK  OF  A 


27 

26 


COMMON  0(  25,2  5) , S ( 25 , 25 ) , I NNODE ( 25 ,25) ,NONOOE(25)  , 
1RH(25,25) ,LH( 25 , 25 ,L INKSl 25 ) , THAT A (25 ,25 ) ,NTOURS, 

1 1 END, NC.TOTO IS, LENGTH (25) , RANGE, ISTOP,?RED(  25,25) 
DIMENSION  X(25) , Y( 25 ) , STAR ( 25 , 25 ) ,NNST R ( 25 , 25 ) , 
1LLSTR(25,25) ,S£CSTR{25,25) ,SNNSTR(25,25) , SL LSTR( 25, 25 ) 
l,CRIT(25) ,OLORH(25,25) ,0LDLH(25,25),Z(25 I 
INTEGER  S,RH,OLORH,OLOLH,EXNOOE,PRED 
REAL  LENGTH 

IPRINT=0 

LIMIT*LINKS(NTOURS) 

DO  26  J>1,KK 
DO  27  1=1, LIMIT 
NL=LH(NTOURS,I) 

NR=RH(NTOURS, I ) 

IN=N0N0DE(J) 

THETA( I ,J)=0(NL,IN)+O(IN,NR)-O(NL,NR) 

CONTINUE 

CONTINUE 


IF  USER  DESIRES 
SET  IPRINT=6 


TO  PRINT  COMPUTED  VALUES  OF  THETA(I,J) 


IFdPRINT  .NE.  5)  GO  TO  50 
DO  112  1=1, LIMIT 

WRITE(6,113)  ( THETA( I,J) ,J=1,KK) 
113  F0RMAT(9F10.2) 

112  CONTINUE 
50  RETURN 
END 

SUBROUTINE  THREEX 


SUBROUTINE  CHECKS  CITIES  IN  A GIVEN  SUBTOUR  TO 
DETERMINE  IF  THE  LENGTH  OF  THAT  SUBTOUR  MAY  BE 
DECREASED  BY  RELOCATING  THE  CITY  IN  THAT  SUBTOJR 

COMMON  0( 25,25) , S ( 25 ,25 ) , I NNODE ( 25 ,25 ) ,NONOOE(25) , 
1RH(25,25) ,LH( 25 , 25 , L INKS! 25 ) , THATA (25 ,2 5 ) ,N TOUR S, 
1IEND,NC,T0TDIS,LENGTH(25) , RANGE, I STOP , PRED ( 25 , 25 ) 
DIMENSION  X(25) , Y( 25 ) , STAR ( 25 , 25 ) ,NNSTR(25  ,25) , 

1LLSTR( 25,25) , SECSTR ( 25 , 25 ) , SNNSTR ( 25 , 25 ) , SL LSTR ( 25, 25 ) 
1,CRIT(25) ,0LDRH(25,25) , OLDLH( 25 , 25 ) , Z( 2 5 ) 

INTEGER  S,RH,aLORH,OLDLH,EXNODe,PRED 
REAL  LENGTH 

55  DO  60  M=1,NT0JRS 

IF(LINKS(M)  .LT.  ^)  GO  TO  ^^8 
LIMIT=LINKS(M)-1 
DO  41  1=1, LIMIT 
NL=LH(M,I  ) 

EXNODE=RH(M,I ) 

NR=S(M,EXNOOE) 

TSTVAL=0(NL,£XNQDE)+D(EXNODE,NR)-D(NL,NR) 

INLIM»I+1 
NEWLIM=LINKS( M) 

DO  12  K=1,NEWLIM 

IF(K.EQ.I.OR.K.EQ.INLIM)  GO  TO  12 
IL  = LH(M,K) 

IR»RH(M,K) 

COMVAL*D( IL,EXN00E ) +0 ( EXNO DE , I R )-D( I L , IR ) 

IF(COMVAL  .LT.  TSTVAL)  GO  TO  45 
12  CONTINUE 


h3 


41  CONTINUE 

48  IF(M  .EQ.  NTOURS)  GO  TO  46 
60  CONTINUE 
GO  TO  46 

C CREATE  THE  IMPROVED  TOJR 

45  S(M,IU»EXNOOE 

S(MtEXNOOE)»  IR 
S (MtNU  -NR 
CALL  ENOPTS 
CALL  TOROIS 
GO  TO  55 

46  RETURN 
END 

SUBROUTINE  TOJRNU 

C SUBROUTINE  CREATES  A NEW  SUBTOUR  WHEN  IT  IS  NOT 

C FEASIBLE  TO  COMPLETE  A VISIT  TO  EACH  CITY  WITH  THE 

C PRESENT  NUMBER  OF  SUBTOURS 


COMMON  0(25.25) « S( 25*25) tINNOOE (25*25) tNQNOOE (25) * 
1RH(25,25) ,LH(25,25*LINKS(25) , THATA ( 25 ,25 ) .NTOURS , 
II£NO,NC,TOTOIS,LENGTH(25),RANGE, IST0P,PRE0(25,25) 
DIMENSION  X(25) ,Y(25), STAR! 25,25) ,NNSTR( 25,25) , 
ILLSTR(25,25) , S£CSTR( 25 *25 ) , SNNSTR ( 25 , 25 ) , SL LSTR ( 25 , 25 ) 
l*CRIT(25 ),0L0RH(25,25) *0LDLH(25,25)i2(25 ) 


Integer  ^.rh, 
REAL  length 


* ui.uL.n  i I *4.1 

OLORHfOLDLH, EXNOOE.PREO 


AMIN)  GO  TO  63 


AMIN=99999. 

LIMIT»LINKS(NTOURS) 

NEWTOR=NTOURS*1 
00  62  J»I,KK 

N0D»N0N00E( J) 

IF(0(l*NOO)  .LT. 

GO  TO  62 
63  AMIN-0(1*N00) 

IN=NOO 
62  CONTINUE 

S(NEWTOR,l)=IN 

S(NEWTOR,IN)=l 

LINKS(NEWTOR) =2 

NTOURS*NEWTOR 

RETURN 

END 

SUBROUTINE  CRSEXC 

C SUBROUTINE  CHECKS  CITIES  IN  A GIVEN  SUBTQUR  TO 

C DETERMINE  IF  TOTAL  DISTANCE  TRAVELED  MAY  BE  DECREASED 

C BY  PLACING  THEM  IN  A DIFFERENT  SUBTOJR 


COMMON  0(25,25) ,S (25,25 ), INNOOE (25, 25), NONODE (25) , 
1RH(25,25),LH(25,25,LINKS(25) , THATA ( 25 ,25 ) .NTOURS* 

1 1 END, NC, f QTOI S.LENGTH(25) .RANGE, I STOP, PREO(  25,25) 
DIMENSION  X( 25) ,T(25),STAR( 25,25) .NNSTR( 25,25), 
1LLSTR(25,25) ,SECSTR(25,25),SNNSTR(25,25) , SL LSTR ( 25 , 25 ) 
i,CRIT(25) ,0LDRH(25,25) * OLDLH ( 25 , 25 ) , Z( 25 ) 

INTEGER  S,RH,OLORH,OLDLH, EXNOOE.PREO 
REAL  LENGTH 

36  DO  80  I«l. NTOURS 
NOO»S( I ,1) 

PRED(I,N00)»1 
81  INT*S(I.NOO) 

PREDd  , 1NT)=N00 
NOD* I NT 


ih 


IFdNT  .ME.  1)  GO  TO  81 
80  CONTINUE 

00  82  I=1,NT0JRS 
LIMIT*UNKS(  I) 

DO  83  J»1,NT0JRS 

IF( I .EQ.  J)  GO  TD  83 
NJNBER=LIN<S( J )-l 
IF(NUMBER  .:0.  1)  GO  TO  82 
DO  84  K=l, NUMBER 
DO  85  M=l, LIMIT 
NL»LH( I,M) 

NRsRHU.M) 

IN^INNODEI J,K) 

IBFR:^PReO<Jt  IN» 

IAFT=S( Jf IM) 

TMeTA(M,KJ»D(NLf  IN  »+D ( I , NR » -D (NL , NR  ) 
BETA=0(1BFR, I N ) +0 ( IN ♦ I AFT ) -D ( IBFR, I A FT) 
IFITHETAIM.O  .GE.  BETA)  GO  TO  85 
CHECK=THETA(M,K)+LENGTH( I ) 

IFICHECK  .3T.  RANGE)  GO  TO  85 


S( J,IBFR)»IAFT 
S( J,IN)=0 

LINKS! J)*LINKS(J)“1 


INSERT  IN  TOUR  I 


S( I,NL)=IN 
S( It IN)»NR 
LINKS!  DELINKS!  U+1 
CALL  ENOPTS 
CALL  TOVIST 
CALL  TORDIS 
GO  TO  86 
85  CONTINUE 
84  CONTINUE 
83  CONTINUE 
82  CONTINUE 
RETURN 
END 


■i 
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